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We address the influence of tlie molecular orbital geometry and of the molecular alignment with 
respect to the laser-field polarization on laser-induced nonsequential double ionization of diatomic 
molecules for different molecular species, namely N2 and Li2. We focus on the recollision excitation 
with subsequent tunneling ionization (RESI) mechanism, in which the first electron, upon return, 
promotes the second electron to an excited state, from where it subsequently tunnels. We show that 
the electron-momentum distributions exhibit interference maxima and minima due to the electron 
emission at spatially separated centers. We provide generalized analytical expressions for such 
maxima or minima, which take into account s p mixing and the orbital geometry. The patterns 
caused by the two-center interference are sharpest for vanishing alignment angle and get washed out 
as this parameter increases. Apart from that, there exist features due to the geometry of the lowest 
occupied molecular orbital (LUMO), which may be observed for a wide range of alignment angles. 
Such features manifest themselves as the suppression of probability density in specific momentum 
regions due to the shape of the LUMO wavefunction, or as an overall decrease in the RESI yield 
due to the presence of nodal planes. 



A. Introduction 

Strong-field phenomena such as high harmonic gener- 
ation (HHG) or above-threshold ionization (ATI) have 
been used as tools for the attosecond imaging of molec- 
ular orbitals pj], for probing the structural changes in 
molecules with attosecond precision and for studying 
quantum interference effects due do photoelectron or 
high- harmonic emission at spatially separated centers Q ■ 
This has been made possible due to the fact that both 
phenomena are caused by the rescattering or recombina- 
tion of an electron with its parent molecule, which, for 
typical intense lasers, occur within hundreds of attosec- 
onds. The simplest targets for which this interference can 
be studied are diatomic molecules, which can be viewed 
as the microscopic counterpart of a double-slit experi- 
ment ^j. 

Potentially, laser-induced nonsequential double ioniza- 
tion (NSDI) can also be employed for probing molecular 
orbitals since laser-induced recollision plays an impor- 
tant role in this case. In NDSI the returning electron 
rescatters inelastically with its parent ion, or molecule, 
giving part of its kinetic energy to a second electron. 
This electron can be released in the continuum either 
through electron-impact ionization j^-[T2| or recollision 
excitation with subsequent tunneling ionization (RESI) 
[13I . The former recollision mechanism happens when 
the first electron, upon return, gives enough energy to 
the second electron of the target so that it can overcome 
the second ionization potential and reach the continuum. 
The latter recollision mechanism happens when the first 
electron, upon return, gives just enough energy to the 
second electron so that it can be promoted to an excited 
bound state, from where it subsequently tunnels. 

In principle, NSDI exhibits several advantages, with 
regard to ATI or HHG. First, it allows one to ex- 



tract more dynamic information about the system, as 
the type of electron-electron interaction can be identi- 
fied in the electron-momentum distributions 0, [Tl| . 
Furthermore, different rescattering mechanisms, such as 
electron-impact ionization or RESI, populate different re- 
gions in momentum space and hence can also be traced 
back from such distributions flBj. Apart from that, 
events happening at different half cycles of the driv- 
ing field can be mapped into different momentum re- 
gions. Concrete examples are NSDI with few-cycle pulses 
[l6|, which lead to asymmetric electron- momentum dis- 
tributions, and individual processes in NSDI of diatomic 
molecules ■ Finally, electron-electron correlation is at 
the essence of this phenomenon and cannot be ignored. 
In contrast, for high-order harmonic generation one may, 
to first approximation, only consider a single active elec- 
tron and the highest occupied molecular orbital. In fact, 
only very recently have multiple orbitals and electron- 
electron correlation been incorporated in the modeling 
of molecular high-order harmonic generation (T8l - [2^ . 

For the above-mentioned reasons, NSDI of molecules 
is being increasingly investigated since the past few 
years. In fact, there has been experimental evidence 
that the orbital symmetry [23j and the alignment an- 
gle 24] affect the shapes of the electron momentum dis- 
tributions. Since then, many theoretical studies have 
also been performed for molecules, involving, for in- 
stance, classical trajectory methods [2^ . the numeri- 
cal solution of the time-dependent Schrodinger equa- 
tion in reduced-dimensionality models [26j . and semi- 
analytical approaches based on the strong-field approxi- 
mation [l3, [23, HI] . Semi-analytical models for NSDI in 
molecules, however, focus on the electron-impact ioniza- 
tion rescattering mechanism. 

For instance, in our previous paper [2^ we addressed 
the infiuence of the orbital symmetry and the molecular 
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alignment with respect to the laser-field polarization on 
NSDI of diatomic molecules for the electron-impact ion- 
ization mechanism. We showed that the electron momen- 
tum distribution exhibit interference maxima and min- 
ima due to electron emission at spatially separated cen- 
ters. Such fringes were positioned at +P2|| = const., 
i.e., parallel to the anti diagonal of the plane spanned by 
the electron momentum components p„|| n = 1, 2 parallel 
to the laser-field polarization. They were sharpest if the 
molecule was aligned along the direction of the field, i.e., 
for vanishing alignment angle. As this angle increased, 
the fringes got increasingly blurred until they were com- 
pletely washed out for perpendicular alignment. 

Apart from that, recently, several studies have found 
that the core dynamics, in particular excitation, is impor- 
tant for high-harmonic generation in molecules (20l.l2^. 
and in particular for attosecond imaging of matter. We 
expect this also to be the case for nonsequential double 
ionization. For that reason, in the past few years, we 
have focused on the RESI mechanism. We have shown 
that the shape of the electron momentum distributions 
depends very strongly on the initial and excited bound 
states of the second electron [H, [s^l, in fact far more 
critically than for electron-impact ionization Q. If this 
is the case already for single atoms, one expects this de- 
pendence to be even more critical for molecules. 

For RESI, we expect the electron momentum distri- 
butions to be affected very strongly by the geometry of 
the bound-state wavefunctions, not only because the ex- 
citation process strongly depends on them, but also due 
to the fact that the second electron is reaching the con- 
tinuum by tunneling. It is by now well known that this 
ionization mechanism is strongly influenced by the pres- 
ence of nodal planes or the directionality of a particular 
molecular orbital. For instance, for HHG the nodal plane 
of a TT state suppresses tunnel ionization when it coincides 
with the polarization axis (see, e.g., [isl. [20l. [2^ Isil - lssj ) . 

In the present paper, we perform a systematic analy- 
sis of quantum-interference effect in NSDI of diatomic 
molecules considering the RESI mechanism. We con- 
struct a semi-analytical model, based on the strong-field 
approximation (SFA), in which an electron tunnels from 
the HOMO of a neutral molecule and rescatters with the 
HOMO of its singly ionized counterpart. Thereby, we 
assume that the second electron is excited to the lowest 
unoccupied molecular orbital (LUMO). We investigate 
the influence of such orbitals and of the alignment angle 
on the NSDI electron momentum distributions. Specif- 
ically we choose species for which these orbitals have 
different geometries and parities. Furthermore, we ad- 
dress the question of whether well-defined interference 
patterns such as those observed in ATI or HHG compu- 
tations may also be obtained for NSDI in the context of 
the RESI mechanism, and, if so, under which conditions. 
These are complementary studies to those performed in 
our recent work on RESI [29, 30], where we show that, for 
single atoms, the shapes of the electron momentum dis- 
tributions carry information about the bound state from 



which the second electron leaves and the state to which 
it is excited. 

This paper is organized as follows. In Sec. U we dis- 
cuss the expression for the RESI transition amplitude, 
including its general expression (Sec. II Ap . the saddle- 
point equations obtained from it (Sec. II B[) and the spe- 
cific prefactors for a diatomic molecule using Gaussian 
orbital basis sets (Sec. II CI) . At the end of this section, 
fSec. lID]) we derive a general two-center interference con- 
dition for the RESI mechanism. Subsequently, in Sec. [Hi 
we compute electron momentum distributions, with em- 
phasis on the two-center interference fSec. lII At . and the 
influence of different molecular orbitals (Sec. Ill Bp . Fi- 
nally, in Sec. IIIII we state the main conclusions to be 
drawn from this work. 



I. STRONG-FIELD APPROXIMATION 
TRANSITION AMPLITUDE 

A. General expressions 

The SFA transition amplitude describing the RESI 
mechanism reads (for details on the derivation see [SQ]). 

/OO pt pt p 

dt dt' dt" d^k 
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and 



B. Saddle-point analysis 



(P. Wilson l^f))^^ 



(27r)3/2 



X / dV2yio„(r2)cxp[-zp2(<)-r2]7Af^(r2).(5) 



Eq. ([T|) describes the physical process in which, at 
a time t", the first electron tunnels from a bound state 

IV'g^'' > into a Volkov state |k(i') >. Then the released 
electron propagates in the continuum from t" to t' , and 
it is driven towards its parent molecule. Upon return, 
the electron scatters inelastically with the core at t' and, 
through the interaction V12 , promotes the second electron 
from the bound state \'4>g > to the excited state j^/^e >• 
Finally, at a later time t, the second electron, initially in a 
bound excited state >, is released by tunneling ion- 
ization into a Volkov state |p2 (i) >. In the above-stated 
equations, Eng {n — 1,2) are the ionization potentials of 
the ground state, Ene {n — 1,2) denote the absolute val- 
ues of the excited-state energies and the potentials Vb(ri) 
and Mon(i"2) correspond to the neutral molecule and the 
singly ionized molecular species, respectively. Here, the 
final electron momenta are described by p„(n = 1, 2). All 
the information about the binding potentials viewed by 
the first and second electrons and the electron-electron 
interaction are embedded in the form factors ([3]) , ([5]) and 
(HI) respectively. Assuming that the electron-electron in- 
teraction depends only on the difference between the two 
electron coordinates, i.e., if Vi2(ri_r2) = Vi2(ri_r2), Eq. 
Q may be rewritten as 
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(27r)3/2 

rf3^^g-»(pi-k).r.[^(2)(j.^)]*^(2)(j.^)^(g) 



with 

VMpi - k) = J d^re-^^P^-'^y^VMr) (7) 

and r = ri_r2. 

Within the framework of the SEA these prefactors 
are gauge dependent. Specifically, in the length gauge 
Pn (r) = Pn + A(r) and k(T) = k + A(r)(T_= t',t"), 
while in the velocity gauge p„ (r) = p„ and k(T) = k. 
This is due to the fact that the gauge transformation can- 
cels out with the minimal coupling in the latter case. In 
practice, however, for the specific situation addressed in 
this work, both gauges lead to very similar results. This 
happens as the above-stated phase differences will cancel 
out in Vpie.kg, and, in Vp^e, ■A-(t) ~ for the parameter 
range of interest (for more details see [1^). In the fol- 
lowing, unless strictly necessary, we will drop the time 
dependence in p„ (r) . 



Subsequently, the transition amplitude ([T]) is solved 
employing saddle-point methods. Eor that purpose, 
one must find the coordinates {ts,t'g,t'^,'k.s) for which 
5(p„, k, t, t', t") is stationary, i.e., for which the con- 
ditions dtS{pn,k,t,t',t") = dt'S{pn,k,t,t',t") = 

df"S{pn,k,t,t',t") = and dwS{pn,k,t,t' ,t") = are 
satisfied. This leads to the equations 



[k-t- A(t")] = -2E_ 



l9' 



(8) 



(9) 



[pi + A{t')f = [k + A(t')]' - 2{E2g - E2e). (10) 



and 



[P2 + A(t)] 



-2Eo 



(11) 



which, as discussed below, provide additional physical 
insight into the problem. 

The saddle-point Eq. ([SJ gives the conservation of en- 
ergy when the first electron tunnel ionized at a time t" . 
Eq. ([9]) constraints the intermediate momentum k of the 
first electron and it makes sure the electron returns to the 
side of its release, which lies at the geometrical center of 
molecule. Eq. (|10|) expresses the conservation of energy 
at a time t' , when the first electron rescatters inelastically 
with its parent ion, and gives part of its kinetic energy 
Eret{t') = [k + A{t')]^ /2 to the core to excite the second 
electron from a state with energy E2g to a state with en- 
ergy i?2e- Immediately after rescatering the first electron 
reaches the detector with momentum pi. Einally, Eq. 
([TT|) describes the fact that the second electron tunnels 
at time t from an excited excited state E2e and reaches 
the detector with momentum p2. As a consequence of 
the fact that tunneling has no classical counterpart, these 
equations possess no real solutions (for more details see 



C. Molecular prefactors 

In this work, we consider that all molecular orbitals 
are frozen apart from the HOMO and the LUMO. We 
also assume frozen nuclei and a linear combination of 
atomic orbitals (LCAO) to construct approximate wave 
functions for the active orbitals. This implies that the 
molecular bound-state wave function for each electron 
reads 

V(")(r„) = ^c„[0(")(r„+R/2)+(-l)'''+^°e)(r„-R/2) 



(12) 
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where R and la denote the internuclear separation and 
the orbital quantum numbers, respectively. The index 
n — 1,2 refers to the electron in question. The index 
Act = applies to gerade symmetry and Aq, = 1 to unger- 
ade symmetry. The binding potential of this molecule, 
as seen by each electron, is given by 



V^Vn) = V^(r„ - R/2) + V^(r„ + R/2) 



(13) 



where the subscript x = or ion refers either to the 
neutral molecule or to its ionic counterpart, respectively, 
and Vjf(r„) = Zcff/r„ is the potential at each center in 
the molecule. Thereby, ZeS is the effective core charge 
as seen by each of the two active electrons. 



In this paper the wave function 
by a Gaussian basis set. 



is approximated 



I la la la 



exp[-C,r 



(14) 



The coefficients bj and and the exponents Q can be 
extracted either from existing literature or from quantum 
chemistry codes. We compute these coefficients using 
GAMESS-UK [13. In our basis set, we took only s and 
p states. This means that, in all the expressions that 
follow, la and Ijs are either or 1. 

The above-stated assumptions lead to the form factors 
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and 
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^2, 

(17) 



l2= / rfV2Vo(r2)e-^P-'-^0i2)(j.^)^ (^8) 



where 



In general, the form factor (jS)) does not affect the shape 
of the electron-momentum distributions. This is particu- 
larly true when the first electron tunnels from an orbital 
with no nodal planes, such as a Cg orbital (2^. However, 
one has to be careful when the electron tunnels from any 
orbital with at least one nodal plane, such as a tt or- 
bital, as this would lead to a suppression of ionization 
for specific alignment angles. 




FIG. 1; Schematic representation of the molecule and laser 
field frames of reference, represented by the black and red sets 
of axis X, y, z and x' , y' , z' respectively. The two centers of the 
molecule are apart by R along the z axis of the molecule, and 
their positions are indicated by the blue circles in the figure. 
The field A(t) is polarized along the z' axis, and 6 shows the 
alignment angle of the molecule with respect to the laser field. 



In the following, we will write the above-stated equa- 
tions as functions of the electron-momentum components 
Pn\\ and Pn± parallel and perpendicular to the laser- field 
polarization. Physically, we are investigating a diatomic 
molecule whose main axis is rotated of an angle 6 with 
respect to the direction of the laser-field polarization. 
Hence, we are dealing with two frames of reference, i.e., 
the molecular frame of reference and the laser field frame 
of reference. The electron momenta in terms of their par- 
allel and perpendicular components with regard to the 
laser-field polarization read 



Pn =Pn\\ez' + Pn± COS (pea;' + Pn± SilKfiBy', 



(19) 



where we assumed that the laser field is polarized along 
the z' axis, the coordinates x' and y' define the plane 
perpendicular to the laser-field polarization and cp is the 
azimuthal angle. In order, however, to compute the 
momentum-space wavefunctions for this molecule, we 
need the momentum coordinates in the frame of refer- 
ence of the molecule. The molecular coordinates x, y and 
z can be obtained by a coordinate rotation around the x 
axis. In this case, the momenta of the electrons in terms 
of parallel and perpendicular components in this latter 
frame of reference will be 



Pn = (p„|| cosS +p„^ sin6'sin(/5)ez +p„^ cos(pea; 

+ (Pni cos 6* sin (/3 — Pn\\ sm6)ey. (20) 

This implies that the momentum components Prix^Pny 
and pnz are defined by Eq. ((20)) and that 

p„ • R/2 = (p„|| cos6' +Pni_ sin6'sin(^)i?/2. (21) 
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A schematic representation of both the field and molec- 
ular sets of coordinates is presented in Fig. [T] Be- 
low, we provide the explicit expressions for the integrals 
1,1 (n — 1,2) in the prefactors (IT5|) and ([T7|. for the spe- 
cific types of orbitals employed in this work. 



1. Excitation a a 

If the second electron is excited from a a to a, a orbital, 
both integrals will have the forms 



I^ = 



E 



where 



and 



where 



1 2'°+'^(Cj +0')^/'+'°+'" 
xexp[- ]P^;^^J^ ].TG.,;,) (22) 



1, la+lp=0 
(pi-k)^, la+ll3^l 

2(0+C/)-(Pi-k)2, L + lp^2 



I2 = Y.bf{-^y^G{lp), 
3' 



G(h) = 



(p2./p2)/^^'\;/3 = l 



(23) 
(24) 

(25) 



In Eq. (|25l) . Ir" and Ir" indicate the radial 
integrals 



One should note that, if the electron is excited from 
a TT to a tr orbital, Ii will also have this form. In the 
second pref actor. 



P2 



(28) 



with Ip — 1. Throughout, (pi — h)^ and with x 
x,y,z are defined according to Eq. (^0)) . 



D. Interference Condition 

Here we provide a general interference condition, which 
takes into account the structure of the orbitals. This in- 
cludes s p mixing and the orbital parity. The expressions 
that follow are easily derived if the exponentials in Eqs. 
([T5l) and ([TTl) are expanded in terms of trigonometric 
functions. In this case, the prefactor ([15]) can be written 
as 



Pie, kg 



with 



^tIw^ E E y^r^smKr + (pr-k).R/2], 

(29) 



a p 



S^i — arctan[ 



and 



C± = i±(-i)'=+'/3+>'=+^^ 



(30) 



(31) 



A similar procedure for high-order harmonic generation 
has been adopted in Sli] . Interference minima are present 
if 



r^^'^^Jhihr) exp[-Or']dr, (26) 



where ji„ {ji2f) denotes spherical Bessel functions. 



Ci + (pi - k) • R/2=m^, 



(32) 



where m is an integer. Similarly, interference maxima 
are obtained for 



2. Excitation a — tt 

We also consider that the second electron is excited 
from a a orbital to a tt orbital. In this case, these orbitals 
are degenerate. For that reason, we choose to consider a 
coherent superposition of the and tt^ orbitals carrying 
equal weights. This gives 



J-j' 

H(Pi-k).)'- (Pi-k)^ 

2'"+'''(0+C,')3/^+'°+'^ ^ 4(0+0')^- ^ ^ 



-f (Pi - k) • R/2=(2m + l)7r/2. (33) 

We will focus on the minima given by Eq. ([32]) as they 
are much easier to observe. If this equation is written in 
terms of the electron momentum component (pi — \i)z 
parallel to the molecular axis we find 

[(pi|| — k) cos 6 + pi± s'm 9 sin ip\ R/2 = rmr — ^1. (34) 

The above-stated equation shows that the parallel mo- 
mentum component pi\\ parallel to the laser-field polar- 
ization will lead to well-defined interference fringes ap- 
proximately at 

2(m7r-Ci) 

Pill = — 5 a -^5) 

" RcosO 
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This means that, in the plane Pi\\P2\\, these minima will 
be at pi 1 1 = const., i.e., parallel to the axis. The per- 
pendicular component pi± will mainly cause a blurring 
in such fringes, when the azimuthal angle is integrated 
over. Extreme limits will be found for alignment angle 
9 = 0, with sharp two-center patterns, and 9 = 90°, when 
they get washed out. 

Following the same line of argument. 



with 



and 



47r 



(27r)3/2 



E 



sin[e2+P2-R/2]X2, (36) 



^2= arctan[ 



(37) 



i:>± 1 ± (-1)'"+^" . (38) 

Interference minima are present for Eq. p6l) if 



^2 + P2 • R/2=TO7r 

Likewise, there will be interference fringes for 



P2\\ 



2(m7r-6) 
R cos 9 



(39) 



(40) 



i.e., parallel to the p2\\ axis in the plane spanned by the 
parallel momentum components Pi\\, P2\\- In the veloc- 
ity and the length gauges, P2\\ = P2\\ and P2\\ + A(t), 
respectively. Since, however, A{t) ~ for the electron 
tunneling time, in practice there will be very little dif- 
ference. The perpendicular momentum components will 
lead to a blurring in the fringes. 



II. ELECTRON MOMENTUM DISTRIBUTIONS 

In this section, we compute electron momentum dis- 
tributions, as functions of the momentum components 
{Pi\\tP2\\) parallel to the laser-field polarization. We as- 
sume the external laser field to be a monochromatic wave 
linearly polarized along the axis z' . Explicitly, 



F,{t) = eo sinujiiz 



(41) 



This approximation is reasonable for laser pulses of the 
order of ten cycles or longer Q. These distributions, 
when integrated over the transverse momentum compo- 
nents, read 

FiPi\\,P2\\) = JJ<fpi±d^P2±\MR{pi,p2) (42) 

+ Ml(pi,P2) + pi ^ P2p, 

where Af/?(pi,P2) and Afi(pi,p2) refer to the right and 
left peak in the electron momentum distributions, respec- 
tively, the transition amplitude Mji{pi,p2) is given by 



Eq. ([Ij, and cPpni. = Pn±dpn±dipp^ . For a monochro- 
matic field, we can use the symmetry A(t) = —A{t ± 
T/2), where T = 27r/aj corresponds to a field cycle, in 
order to simplify the computation of the electron mo- 
mentum distributions. This is explained in detail in 
our previous work 29]. We also symmetrize the above- 
stated distributions with respect to the particle exchange 
Pi O p2. To a good approximation, it is sufficient to con- 
sider the incoherent sum in Eq. (|42l) as the interference 
terms between the right and left peaks practically get 
washed out upon the transverse momentum integration 
(see Appendix B in (29|). 

In the following, we will compute electron momentum 
distributions for Li2 and N2. For all cases, we assume 
that the electron-electron interaction is of contact type, 
i.e., V12 = 6{ri — r2). This will avoid a further momentum 
bias in the electron-electron distributions as it leads to 
^12 (pi ~ k) = const, and allow us to investigate the 
influence of the target structure alone. For a long-range 
potential, Vi2(pi — k) would be momentum dependent, 
and hence mask the features we intend to investigate. 



A. Interference effects and s p mixing 

We will commence by investigating whether the inter- 
ference conditions derived in Sec. lIDI hold. For that pur- 
pose, we must have non-negligible tunneling ionization 
for parallel-aligned molecules, as this is the situation for 
which the fringes are expected to be sharpest. Hence, one 
must consider a target for which neither the HOMO nor 
the LUMO exhibits nodal planes along the internuclear 
axis. Therefore, we assume that the first electron tun- 
nels from the HOMO in Li2 and rescatters inelastically 
with Li^, exciting the second electron from its HOMO 
(2(Tg) to its LUMO (2ct„). In order to get a clear pic- 
ture of conditions p2|) and ([39|) . we must investigate the 
corresponding prefactors individually. 

In Fig. [21 we depict the above-mentioned electron- 
momentum distributions for alignment angle 9 = 0°. We 
consider V^^e.kg = const, and focus on the influence of 
V^2e alone. We take either the individual contributions 
of s and p states or the combination of both for 2(t„. For 
clarity, in the upper panels, we also exhibit the distri- 
butions obtained without symmetrizing with respect to 
the momentum exchange and electron start times. For 
all cases, the two-center fringes in Fig. [2] are parallel to 
P2\\ = const., in agreement with the second interference 
condition derived in Sec. II Dl 

For pure s or p states and = 1, which is the case 
for a Uu orbital, this condition can be further simplified. 
It reduces to 



for s states, and 



sin[p2 • R/2] = 0, 



cos[p2 • R/2] = 



(43) 
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FIG. 2: Electron-momentum distributions for NSDI in 
Li2 (bound-state energies Eig — 0.18092040 a.u., E2g — 
0.43944428 a.u. and E2e = 0.12481836 a.u. and equilib- 
rium internuclear distance R = 4.7697 a.u.) considering 
only the RESI mechanism, as functions of the momentum 
components parallel to the laser-field polarization, obtained 
considering Vpje according to Eq. ((Sj and Vpje.kg ~ const. 
We consider zero alignment angle, driving- field intensity / — 
4.6 X lO^^W/cm^ and cj = 0.057 a.u. respectively. Panels (a) 
to (c) display only the contribution from the orbits starting 
in the first half cycle of the field, while in panels (d) to (f) 
the distributions have been symmetrized to account for the 
electron orbits starting in the other half cycle and for elec- 
tron indistinguishability. The left, middle and right panels 
correspond to the contributions of the s, p and all states used 
in the construction of the au LUMO, respectively. The solid, 
dashed and short dashed lines show the position of minima 
due to the two-center interference, node of the wavefunction 
and mixed cases, respectively. The contour plots have been 
normalized to the maximum probability in each panel. 



for p states. This implies that, for the former, we expect 
minima at p2 • R = 2m7r, while for the latter they should 
occur at P2 • R = {2m+ 1)tt. The position of such minima 
can also be determined analytically by considering that 
the second electron tunnels at the peak of the laser field, 
i.e., at t ~ tt/2. The dashed lines in the figure show that 
the position of these minima exhibit a very good agree- 
ment with this simple estimate. Physically, this good 
agreement may be attributed to the fact that the second 
electron tunnels most probably at this time. 

For the s states the two-center interference gives a 
sharp minimum at p2\\ = (Figs.[2](a) and (d)), while for 
the p states these patterns are located near p2\\ = iiWUp 



I 

(Figs. [21(b) and (e)). In the p— state case the distribution 
has another minimum at P2\\ = 0, which comes from the 
fact the p wavefunctions vanish for p„ =0. This causes a 
suppression in the transition amplitude. If the contribu- 
tions of both s and p states are considered, the minima in 
the high-momentum region due to the two-center inter- 



FIG. 3: Electron-momentum distributions for RESI in Li2 as 
functions of the electron momentum components parallel to 
the laser-field polarization considering Vpie,kg ~ const and 
Vp2e according to Eq. ((S)), for alignment angles 6 = 45° (pan- 
els (a) to (c)), and 90° (panels (d) to (f)). The remaining 
parameters are the same as in the previous figures. The solid 
lines show the position of minima due to the node of the one- 
center wavefunction. From left to right, we considered the 
contributions of the s, p and all states used in the construc- 
tion of the LUMO. All panels have been symmetrized with 
regard to the electron orbits and indistinguishability. The 
contour plots have been normalized to the maximum proba- 
bility in each panel. 



ference seen for the p states vanish, but the minimum at 
P211 — survives. This is shown in Figs. [21(c) and (f) for 
unsymmetrized and symmetrized distributions, respec- 
tively. 

One should note, however, that for parallel-aligned 
molecules, both the two-center minimum for the s states 
and the minimum caused by the node in the p states 
occur at the same momentum, i.e., at p2|| = 0. Hence, 
when s p mixing is included both mechanisms contribute 
to the suppression at the axes Pn\\ =0 seen in Figs. [21(c) 
and (f). We will now investigate the behavior of this node 
when the alignment angle is varied. Since for Li2 both the 
LUMO and the HOMO exhibit distinct shapes and sym- 
metries one can expect significant changes in the electron- 
momentum distributions when this angle is modified. 

Hence, in Fig. [31 we consider the same prefactors as in 
the previous case, but alignment angles 9 = 45° and 90°. 
The figure shows that the patterns caused by the electron 
emission at spatially separated centers get washed out for 
such angles. This is due to the momentum components 
perpendicular to the laser-field polarization, and can be 
seen very clearly in Fig. [31(a), where the s contributions 
are displayed for 9 = 45°. Already for this angle the 
interference minima at the axes p„|| — are absent. In 
contrast, the suppression at the axes caused by the fact 
that the p wavefunctions vanish in that momentum re- 
gion is still present. This is shown in Fig. [31(b), in which 
the contributions from the p states are depicted. The 
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FIG. 4: RESI electron-momentum distributions for Li2 con- 
sidering Vpje = const, and V^p^e.kg according to Eq. (|15|l . for 
6 = Q. The field and molecular parameters are the same as in 
the previous figure. The upper panels display only the con- 
tribution from the sets of orbits starting in the first half cycle 
of the laser field. In the lower panels the distributions have 
been symmetrized in order to account for the orbits starting 
in the other half cycle of the field, and for electron indistin- 
guishability. The left, middle and right panels display the 
contributions from s, p and all states composing the HOMO 
and the LUMO, respectively. The dashed line shows the po- 
sition of the two-center interference minimum. The contour 
plots have been normalized to the maximum probability in 
each panel. 



blurring is caused by the fact that, in momentum space, 
these wavefunctions are proportional to G{lp — 1) (see 
Eq. (|25p ). This function contains components of P2 both 
parallel and perpendicular to the laser field polarization, 
and the contributions from the latter tend to wash out 
the minimum. When both s and p contributions are con- 
sidered, there is a strong suppression of the yield near 
the Pn\\ = axis (see Fig. 131(c)). We have verified that 
this is due to the destructive interference between both 
types of contributions in this momentum region. 

For 9 — 90°, only the components p2i_ contribute, and 
the electron momentum distributions are determined by 
the momentum-space integration alone. As a result, they 
refiect the momentum-space constraints for the RESI 
mechanism. These constraints lead to electron momen- 
tum distributions peaked at {pi\\,Pj\^) — (±2-^/L'p, 0), 
with j,j — 1,2 and i ^ j and with widths I^JTTp^ and 
have been explicitly written in [29, 30]. This holds both 
for the s, p and mixed case (Figs. 131(d), (e) and (f), re- 
spectively). 

We will now focus on the interference condition de- 
termined by the excitation prefactor With this ob- 
jective, we will keep V^^e = const, and investigate the 
influence of lpie,kg alone, starting from vanishing align- 
ment angle. Once more, we will study the contributions 
of the s and p states, and the overall distributions. The 



interference condition and also the wavefunctions in the 
excitation prefactor now incorporate the HOMO and the 
LUMO (see Eq. (g])). For Li J, the former and the latter 
are a gerade and an ungerade orbital, so that Aq = and 
A/3 = 1 in Eq. p2|) . For pure s states, la = Ip = and for 
pure p states, la — Ip = 1- This will lead to the simplified 
interference condition 



sin[(pi - k) • R/2] = 



(45) 



for both. Hence, one expects a minimum close to vanish- 
ing parallel momenta in the pure cases. When s p mixing 
is included, however, different angular momenta will also 
be coupled and the general interference condition must 
be considered. 

The electron momentum distributions obtained in this 
way are shown in Fig. 21 for both symmetrized and un- 
symmetrized distributions (upper and lower panels, re- 
spectively). For most distributions in the figure, we do 
not observe a clear suppression of the probability densi- 
ties in any momentum region. This holds both for those 
caused by the two center interference and by the geom- 
etry of the wavefunctions at the ions. We have only ob- 
served a two center minimum if we consider the individ- 
ual contributions of the p states, and do not symmetrize 
the distributions (see Fig. 111(b)). This is due to the 
fact that, for the parameters considered in this work, the 
two-center minimum according to condition psp lies at or 
beyond the boundary of the momentum region for which 
rescattering of the first electron has a classical counter- 
part. The center of this region is roughly at ~ 
and its extension is determined by the difference between 
the maximal electron kinetic energy upon return and the 
excitation energy — i?2e , as discussed in our previous 
article [H]. 

Apart from that, s p mixing will lead to a blurring of 
this minimum, as it couples states with different angular 
momenta. Symmetrization introduces other events, ei- 
ther due to the electron indistinguishability or displaced 
by half a cycle, and obfuscates this minimum further, as 
shown in the lower panels of the figure. 

If the alignment angle is varied, incorporating only the 
excitation prefactor Vp-^e,kg will lead to ring-shaped dis- 
tributions, regardless of whether only p, s or all basis 
states employed in the construction of the HOMO and 
LUMO are taken. This is expected as, apart from the 
above-mentioned s p mixing, which blur wavefunction- 
specific features, there will now be transverse momentum 
components in the prefactor Ipie.kg which will wash out 
two-center interference patterns. We have verified that 
this is indeed the case (not shown). 



B. Molecular orbital signature 

In this section, we make an assessment of how the ge- 
ometry of the HOMO and the LUMO affect the RESI 
electron momentum distributions. With this objective, 
we incorporate both prefactors Vp^e and Vp^e,kg and vary 
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the alignment angle. In order to discuss the influence of 
nodal planes, we are also providing the overall yield ob- 
tained in our computation in the two figures that follow 
(see color maps on the right-hand side of each panel). 



From other strong-field phenomena, it is well-known that 
the presence of nodal planes may suppress the yield con- 
siderably [ll,|3li. 




-4 -2 




FIG. 5: Electron-momentum distributions for Li2 as functions 
of the parallel momenta (pi||,P2||) considering all prefactors, 
for different alignment angles. Panel (a), (b) and (c) corre- 
spond the alignment angle 6= 0,45 and 90 degrees, respec- 
tively. The field and molecular parameters are the same as in 
the previous figures. 



We will commence by having a closer look at Li2. Such 
results are displayed in Fig. [Sj The main conclusion to 
be drawn from the figure is that the prefactor V^^e plays 
the dominant role in determining the shapes of the elec- 
tron momentum distributions. This can be observed by 
a direct comparison of Fig. [SJ(a) with Fig. [21(f), for van- 
ishing alignment angle. The distributions in both figures 
exhibit similar shapes and minima at the axes = 0, 
and are very different from those obtained if only the 
recollision-excitation prefactor is included (see Fig.Hff)). 
The main effect of the excitation prefactor Vpje.kg is to 
alter the widths of the distributions. This situation per- 
sists for larger angles, such as 9 = 45° and 6 = 90°, as a 
comparison of Figs. [5](b) and (c), with Fig. [31(c) and (f) 
shows. For 9 = 45°, there is a suppression of the yield 
near the axes p^w = 0, while for 9 = 90° the interference 
patterns are washed out. 

Another interesting feature is that the overall yield de- 
creases with the alignment angle between the molecular 
axis and the field. This is due to the fact that the LUMO, 
from which the second electron tunnels, is a cr orbital. 
Spatially, cr orbitals are localized along the internuclear 
axis, and do not exhibit nodal planes for vanishing align- 
ment angle. This implies that tunneling ionization is 
favored when the LUMO is parallel to the laser field, and 
decreases when the difference between the orientation be- 
tween the field and the LUMO increases. 

A legitimate question is, however, how the shape of the 
molecular orbital to which the second electron is excited 
is imprinted on the electron momentum distribution, if 
there are nodal planes parallel or perpendicular to the 



molecular axis. For that reason, we now present electron 
momentum distributions under the assumption that the 
second electron is excited to a tt^ orbital. Specifically, 
we choose N2 and its singly ionized counterpart, i.e., 
as the molecular species in our RESI computation. The 
first electron will be ripped off from the HOMO, which 
is a 3<7g orbital. However, upon return, it will excite the 
second electron to the LUMO, which is a Itt^ orbital. 
A iTTg orbital exhibits two nodal planes, which will be 
oriented along the laser-field polarization for parallel and 
perpendicular-aligned molecules, i.e., at alignment angles 
6* = 0° and 6* = 90°. This orbital also exhibits lobes at 
angles 9 — {2n + 1)ti/A with regard to the internuclear 
axis. 

The results obtained for this molecular species are ex- 
hibited in Fig. [6l As an overall pattern, we observe that 
the NSDI signal no longer decreases monotonically with 
increasing alignment angle. In fact, the signal increases 
for alignment angle Q < 9 < 45°, is strongest for 9 — 45°, 
and decreases once more for 45° < 9 < 90°. This may 
be easily understood as a consequence of the geometry 
of the iTTg orbital. For 9 = Q° (Fig. [HJa)), the external 
field is parallel to one of the nodal planes. Hence, tunnel 
ionization is strongly suppressed. This refiects itself in 
the overall yield. As the alignment angle increases, the 
field-polarization direction gets further and further away 
from the direction of this nodal plane, and the yield in- 
creases until 9 = 45° (Fig. [HJb)). For this angle, the 
field is parallel to one of the lobes of the liTg orbital, so 
that tunnel ionization of the second electron is enhanced. 
As the alignment angle is further increased, the direction 
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FIG. 6: Electron-momentum distributions for N2 (bound- 
state energies E\g — 0.63486 a.u., E^g — 1.12657 a.u., and 



E2. 
R- 

intensity / = 
momenta (pi 
angles ^ = 0. 
spectively). 



= 0.26871290 a.u. and equilibrium internuclear distance 
2.11 a.u.) in a linearly polarized monochromatic field of 
1.25 X 10^"'W/cm^ as functions of the parallel 
1P211) considering all prefactors, for alignment 
45 and 90 degrees (panels (a), (b) and (c), re- 
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of the field approaches the nodal plane at 6* = 7r/2 and 
ionization is further suppressed (Fig. IHJc)). 

Apart from the above-mentioned behavior, we also ob- 
serve a suppression along the axes pn\\ = 0, regardless 



of the alignment angle. This is due to the fact that, in 
position space, tt orbitals vanish at the origin of the co- 
ordinate system. Consequently, their Fourier transform 
vanish for p„|| =0. 



III. CONCLUSIONS 

The results presented in the previous sections illustrate 
the potential of laser-induced nonsequential double ion- 
ization for the attosecond imaging of molecules. This is 
particularly true if the recollision-excitation with subse- 
quent tunneling ionization (RESI) pathway is dominant. 
The computations in this work show that the shapes of 
the RESI electron momentum distributions depend in a 
dramatic fashion on the geometry of the state to which 
the second electron has been excited by the first electron, 
and from which it tunnels. The state in which the sec- 
ond electron is initially bound, i.e., the highest occupied 
molecular orbital (HOMO) of the singly ionized species 
plays only a secondary role. 

Thereby, two main issues are important in determining 
the shapes of the electron momentum distribution: the 
quantum interference caused by the interference due to 
the photoelectron emission at spatially separated center, 
and the geometry of the orbital from which the second 
electron tunnels. 

In order to investigate the first issue, generalized inter- 
ference conditions for the first and second electron that 



take into account s p mixing along the lines of 3lJ have 



been derived, and led to fringes parallel to the p„ 







n = 1 , 2 axes in the plane spanned by the electron mo- 
mentum components parallel to the laser-field polariza- 
tion. These fringes agreed well with analytic estimates, 
but were washed out for relatively small alignment an- 
gles. 

In contrast, the features caused by the orbital geom- 
etry, such as suppression of the probability density near 
Pn\\ — observed for p states were present over a wide 
range of alignment angles. Furthermore, the presence 
or absence of nodal planes manifests itself as the sup- 
pression, or enhancement, of the overall yield with re- 
gard to the alignment angle. We have discussed the dif- 
ferences and similarities between (t„ and tt^ orbitals in 
this context, exemplified by the LUMOs of N2 and Li2. 
These results agree with those reported in the literature 
for phenomena such as high-order harmon ic g eneration 
I20I [22! and above-threshold ionization [2l| . 
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